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Abstract. We probe the cosmic acceleration by using the recently released SNLS3 sample 
of 472 type la supernovae. Combining this type la supernovae dataset with the cosmic 
microwave background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7- 
yr observations, the baryon acoustic oscillation results from the Sloan Digital Sky Survey 
data release 7, and the Hubble constant measurement from the Wide Field Camera 3 on the 
Hubble Space Telescope, we measure the dark energy equation of state w and the deceleration 
parameter q as functions of redshift by using the Chevallier-Polarski-Linder parametrization. 
Our result is consistent with a cosmological constant at la confidence level, without evidence 
for the recent slowing down of the cosmic acceleration. Furthermore, we consider three binned 
parametrizations (w is piecewise constant in redshift z) based on different binning methods. 
The similar results are obtained, i.e., the ACDM model is still nicely compatible with current 
observations. 
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1 Introduction 

The cosmic acceleration has become one of the most important problems in modern cosmology 
since its discovery in 1998 [1]. This implies that around two thirds of components in our 
universe is composed of a mysterious dark energy (DE) [2]. Although numerous theoretical 
DE models have been proposed in the past decade [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] , the nature 
of DE still remains a mystery. Nowadays an important question concerning the cosmic 
acceleration is to determine whether DE is consistent with a cosmological constant [14] [15]. 

A powerful probe of DE is Type la supernovae (SNIa), which can be used as cosmological 
standard candles to directly measure the cosmic expansion. Recently a high-quality joint 
sample of 472 supernovae (SN), the SNLS3 SNIa dataset [16], was released. This SNIa 
sample includes 242 SN at 0.08 < z < 1.06 from the Supernova Legacy Survey (SNLS) 3-yr 
observations [17], 123 SN at low redshifts [18] [19], 93 SN at intermediate redshifts from the 
Sloan Digital Sky Survey (SDSS)-II SN search [20], and 14 SN at z > 0.8 from Hubble Space 
Telescope (HST) [21]. This SNIa sample has been used to study the evolution of DE equation 
of state (EOS) w by Sullivan et al. [22]. However, some other important quantities, such as 
the DE density p^e [23] and the deceleration parameter q [24] , have not been studied in their 
work. 

In this work, we explore the cosmological consequences of the SNLS3 SNIa dataset. In 
addition to the DE EOS w, we also study the evolution of deceleration parameter q. To 
perform a comprehensive analysis including multiple observational techniques, we combine 
the cosmic microwave background (CMB) anisotropy data from the Wilkinson Microwave 
Anisotropy Probe 7-yr (WMAP7) observations [25], the baryon acoustic oscillation (BAO) 
results from the SDSS Data Release 7 (DR7) [26], and the Hubble constant measurement 
from the Wide Field Camera 3 (WFC3) on the HST [27]. 

In 2009, by studying the popular Chevallier-Polarski-Linder (CPL) parametrization 
[28], Shafiello et al. [29] argued that the Constitution SNIa sample [19] appears to support 
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a dynamical DE. Moreover, their analysis implies a possibility: the cosmic acceleration has 
already peaked, and we are currently witnessing its slowing down. It is interesting to test 
whether this strange result is caused by the Constitution SNIa dataset itself. So in this work 
we revisit the CPL parametrization by using the SNLS3 SNIa dataset. 

Moreover, to further explore the cosmic acceleration, we also take into account the 
binned parametrization. This parametrization was firstly proposed by Huterer and Starkman 
[30] based on the principal component analysis (PCA) [30] [31]. The basic idea is to divide the 
redshift range into different bins and setting w as piecewise constant in redshift z. It should 
be mentioned that there are different methods for the optimal choice of the redshift bins 
in the literatures. In [32], Wang argued that one should choose a constant Az for redshift 
slices. This is because that the observables like H and 1/Da (length scales extracted from 
data analysis) are assumed to be constant in each redshift slice in the galaxy redshift survey. 
In [33], Riess et al. proposed another binning method, which has also drawn a lot of attention 
[34]. In this method, the number of SNIa in each bin times the width of each bin is a constant 
(i.e. nAz = const). In addition, in [35], we presented a new binned parametrization method. 
Instead of choosing the discontinuity points Zi by hand, one can treat Zi as models parameters 
and let them run freely. Since all these choices are reasonable, in this work we perform a 
comprehensive analysis and use all the above methods. 

This paper is organized as follows: In Sec. 2, we describe the phenomenological models 
considered here and the method of data analysis. In Sec. 3, we introduce the observational 
data and describe how they are included in our analysis. In Sec. 4 we present the results 
obtained in this paper. In the end, we give a short summary in Sec. 5. In this work, we 
assume today's scale factor ao = 1, so the redshift z = a~ l — 1; the subscript "0" always 
indicates the present value of the corresponding quantity, and the unit with c = K = 1 is 
used. 



2 Models and Methodology 

Standard candles impose constraints on cosmological parameters through a comparison be- 
tween the luminosity distance from observations and that from theoretical models. In a 
spatially flat Friedmann-Robertson- Walker (FRW) universe (the assumption of flatness is 
motivated by the inflation scenario), the luminosity distance di is given by 

, , , l + z f z dz' . . 

with 

e{z) = h{z)/h = [n r (i + zf + n m (i + zf + (i - n r - n m )f{z)] 1/2 . (2.2) 

Here H{z) is the Hubble parameter, Hq is the Hubble constant, Jl m is the present fractional 
matter density, and $l r is the present fractional radiation density, given by [25], 

n r = 7 (1 + 0.2271iV e// ), fi 7 = 2.469 x 10" 5 /i" 2 , N eff = 3.04, (2.3) 

where S7 7 is the present fractional photon density, h is the reduced Hubble parameter, 
and N e ff is the effective number of neutrino species. The DE density function f(z) = 
Pde(z) I Pde(fy is a key function, because a DE parametrization scheme enters in f(z). 
First, we consider the CPL model, which assumes 

w(z) = w + w a —— , (2.4) 
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where wo and w a are constants. The corresponding f(z) is 

f(z) = (1 + zf^o+^a) exp . (2.5) 

Then, we consider the binned parametrization. For this case, f(z) takes the form 

n-l 

f{z n -l < z < z n ) = (1 + z fV+^) J](l + 2 4 ) 3(mi -^ +l) , (2.6) 

2=1 

where W{ is the EOS in the ith redshift bin defined by an upper boundary at Zj. In this work, 
we consider the case of 3 bins, i.e., n = 3. 

As mentioned above, we consider three binned parametrizations based on three binning 
methods. For the first binning method where Az = const, we choose 

Zl = 0.5, z 2 = 1.0. (2.7) 

We will call it "const Az" model hereafter. 

For the second binning method where nAz = const, we choose 

Zl = 0.3, z 2 = 0.73. (2.8) 

In this way we have nAz ~ 65 at the SNIa reshift region < z < 1.4. We will call it "const 
nAz" model hereafter. 

Finally, in the last binning method, we perform a best-fit analysis, and find the following 
choice 

Zl = 0.2, z 2 = 1.31, (2.9) 

can yield a minimal x 2 - I n our analysis, we set the conditions z\ > 0.1, z 2 > z\ + 0.1 and 
z 2 < 1.4. We will call it "free Az" model hereafter. 

In this work, we adopt the x 2 statistic to estimate the model parameters. For a physical 
quantity £ with experimentally measured value (, bs, standard deviation and theoretically 
predicted value the x 2 takes the form 

X| = fc*'"^) 2 . (2.10) 



The total x 2 is the sum of all x|s, i-e. 



\ 2 



£4 (2-11) 



One can determine the best-fit model parameters by minimizing the total x 2 ■ Moreover, by 
calculating A\ 2 = X 2 ~ Xmin' one can determine the la and the 2a confidence level (CL) 
ranges of a specific model. Statistically, for models with different n p (denoting the number 
of free model parameters), the la and 2a CL correspond to different Ax 2 - In Table 1, we 
list the relationship between n p and Ax 2 from n p = 1 to n p = 9. 

In this work, we determining the best-fit parameters and the la and 2a CL ranges by 
using the Monte Carlo Markov chain (MCMC) technique. We modify the publicly available 
CosmoMC package [36] and generate O(10 6 ) samples for each set of results presented in this 
paper. We also verify the reliability and accuracy of the code by using the Mathematica 
program [37]. 
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Table 1. Relationship between number of free model parameters n p and Ax 2 . 









1 

± 


1 




2 


2.30 


6.18 


3 


3.53 


8.02 


4 


4.72 


9.72 


5 


5.89 


11.31 


6 


7.04 


12.85 


7 


8.18 


14.34 


8 


9.30 


15.79 


9 


10.42 


17.21 



3 Observational data 

In this paper, we use the SNLS3 SNIa sample [16], the CMB anisotropy data from the 
WMAP7 observations [25], the BAO results from the SDSS DR7 [26], and the Hubble con- 
stant measurement from the WFC3 on the HST [27]. In the following, we briefly describe 
how these data are included into the x 2 analysis. 

3.1 The SNIa data 

Here we use the SNLS3 SNIa dataset released in [16]. This combined sample consists of 472 
SN at 0.01 < z < 1.4, including 242 SN over 0.08 < z < 1.06 from SNLS 3-yr observations 
[17], 123 SN at low redshifts [18] [19], 93 SN at intermediate redshifts from the SDSS-II SN 
search [20], and 14 SN at z > 0.8 from HST [21]. The systematic uncertainties of the SNIa 
data were nicely handled [16] . The total data of the SNLS3 sample can be downloaded from 
[38]. 

The x 2 OI the SNIa data is 

X| N = Anl T • Cr 1 • Ant, (3.1) 

where C is a 472 x 472 covariance matrix capturing the statistic and systematic uncertainties 
of the SNIa sample, and Ant = nt# — nt mo d is a vector of model residuals of the SNIa sample. 
Here is the rest-frame peak B band magnitude of the SNIa, and m mo d is the predicted 
magnitude of the SNIa given by the cosmological model and two other quantities (stretch 
and color) describing the light-curve of the particular SNIa. The model magnitude m mot ^ is 
given by 

"imod = 51og 10 £>L(z he i,z cm b) - a{s - 1) + f3C + M. (3.2) 
Here T>l is the Hubble-constant free luminosity distance, which takes the form 

/"^cmb dz 

T>L{ZhB\, 2 ZC mb) = (l + 2hel) / n , (3-3) 

JO ) 

where z cm b and Zhel are the CMB frame and heliocentric redshifts of the SN, s is the stretch 
measure for the SN, and C is the color measure for the SN. a and /3 are nuisance parameters 
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which characterize the stretch-luminosity and color-luminosity relationships, respectively. 
Following [16], we treat a and (3 as free parameters and let them run freely. 

The quantity M in Eq. (3.2) is a nuisance parameter representing some combination 
of the absolute magnitude of a fiducial SNIa and the Hubble constant. In this work, we 
marginalize M following the complicated formula in the Appendix C of [16]. This procedure 
includes the host-galaxy information [39] in the cosmological fits by splitting the samples 
into two parts and allowing the absolute magnitude to be different between these two parts. 

The total covariance matrix C in Eq. (3.2) captures both the statistical and systematic 
uncertainties of the SNIa data. One can decompose it as [16], 

C = D s tat + C s tat + C sys , (3.4) 

where T) sta t is the purely diagonal part of the statistical uncertainties, C s t a t is the off-diagonal 
part of the statistical uncertainties, and C sys is the part capturing the systematic uncertain- 
ties. It should be mentioned that, for different a and /3, these covariance matrices are also 
different. Therefore, in practice one has to reconstruct the covariance matrix C for the cor- 
responding values of a and /3, and calculate its inversion. For simplicity, we do not describe 
these covariance matrices one by one. One can refer to the original paper [16] and the public 
code [38] for more details about the explicit forms of the covariance matrices and the details 
of the calculation of Xsn- 

3.2 The CMB data 

Here we use the "WMAP distance priors" given by the 7-yr WMAP observations [25]. The 
distance priors include the "acoustic scale" I a, the "shift parameter" R, and the redshift 
of the decoupling epoch of photons z*. The acoustic scale I a, which represents the CMB 
multipole corresponding to the location of the acoustic peak, is defined as [25] 

Ia = {1 + z*) — r— . (3.5) 

Here Da (z) is the proper angular diameter distance, given by 

and r s (z) is the comoving sound horizon size, given by 

r (z) = J- [ m+Z) da (3 7) 

VZJo a 2 H{a)^l + (30 fe /4fi 7 )a ' 

where £lb an d 7 are the present baryon and photon density parameters, respectively. In 
this paper, we adopt the best-fit values, ^5 = 0.02253/i~ 2 and f2 7 = 2.469 x 10~ 5 /i~ 2 (for 
T C mb = 2.725 K), given by the 7-yr WMAP observations [25]. The fitting function of was 
proposed by Hu and Sugiyama [40]: 

z* = 1048[1 + 0.00124(ft b /i 2 )- a738 ][l + gi (Q m h 2 ) 92 ], (3.8) 

where 

_ 0.0783(^/i 2 )-°- 238 0.560 
91 ~ l + 39.5(^/i 2 ) - 763 ' 92 ~ l + 21.1(lVi 2 ) 1 - 81 " ( ' 
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In addition, the shift parameter R is defined as [41] 



R(z*) = ^n m H$(l + z*)D A (z*). (3.10) 

This parameter has been widely used to constrain various cosmological models [42]. 
As shown in [25], the \ 2 of the CMB data is 

Xcmb = ( x i bs ~ x f l )(C'CMB)ij( x 'j bS ~ X T^ (3-11) 

where X{ = (Ia,R,z*) is a vector, and (CcAfs)y ^ s ^ ne i nverse covariance matrix. The 7- 
yr WMAP observations [25] had given the maximum likelihood values: Ia(z*) = 302.09, 
R{z*) = 1.725, and = 1091.3. The inverse covariance matrix was also given in [25] 

/ 2.305 29.698 -1.333 
(Ccmb) = 29 - 698 6825.27 -113.180 | . (3.12) 
\ -1.333 -113.180 3.414 

3.3 The BAO data 

Here we use the distance measures from the SDSS DR7 [26]. One effective distance measure 
is the Dy(z), which can be obtained from the spherical average [43] 



D v (z) = 



[l + z) 2 D\{z) 



H(z) 



1/3 



(3.13) 



where Da(z) is the proper angular diameter distance. In this work we use two quantities 
do.2 = r s (zd)/Dy (0.2) and do.35 = r s (^)/-Dy(0.35). The expression of r s is given in Eq.(3.7), 
and Zd denotes the redshift of the drag epoch, whose fitting formula is proposed by Eisenstein 
and Hu [44] 

r 

i + h(n b h 2 ) b2 , (3.14) 



Zd 



i29i(n m h 2 ) - 251 



1 + 0.659(fi m /i 2 ) 



0.828 



where 



61 = 0.313(a m /i 2 )" - 419 [1 + 0.607(O m /i 2 ) - 674 ] , (3.15) 

62 = 0.238(fi m /i 2 ) - 223 . (3.16) 

Following [26], we write the x 2 for the BAO data as, 

Xbao = A Pi( C BAo)ii A Pji ( 3 - 17 ) 

where 

A Pi = pf ata - Pi , pf ta = < ata = 0.1905, pf t& = < a 3* 5 a = 0.1097, (3.18) 

and the inverse covariance matrix takes the form 

! / 30124 -17227 \ 
^bao) ~ I _ 1722 7 86977 ) ' 
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3.4 The Hubble constant data 



The precise measurements of H$ will be helpful to break the degeneracy between it and the 
DE parameters [45]. When combined with the CMB measurement, it can lead to precise 
mesure of the DE EOS w [46]. Recently, using the WFC3 on the HST, Riess et al. obtained 
an accurate determination of the Hubble constant [27] 

H = 73.8 ± 2.4km/s/Mpc, (3.20) 

corresponding to a 3.3% uncertainty. So the x 2 of the Hubble constant data is 

2 //i- 0.738\ 2 

** = (-53a-) ■ (3 ' 21) 

3.5 The total \ 2 

Since the SNIa, CMB, BAO and Hq are effectively independent measurements, we can com- 
bine them by simply adding together the x 2 functions, i.e., 

XaU = XSN + XCMB + XBAO + Xh- ( 3 - 22 ) 



4 Results 



As mentioned above, we study the CPL model and the three binned models (const Az, const 
nAz and free Az) in this work. A brief summary of these models are shown in Table 2, where 
the models, model parameters (together with their best-fit values and la uncertainties), and 
Xmin s are gi ven - The nuisance parameters h, a and j3 used in the analysis are actually not 
model parameters with significant meanings, so we do not list them in the table. 

To make a comparison, we also listed the value of x 2 min for the ACDM model. The result 
shows that the inclusion of extra parameters (wo, w a , w\, W2) does not lead to a remarkable 
reduction of the Xmin- This implies that the ACDM model can provide a nice fit to the current 
data, and the consideration of more complex model is not preferred. In fact, the ACDM 
model usually has the best performance under the model-independent consideration [47], 
based on the "information criteria" [48]. In this work, since our main purpose is to explore 
the cosmological consequences of the SNLS3 sample and to probe the cosmic acceleration, 
we will not discuss the topic of the model comparisons. 



Table 2. Summary of the models considered 



Model 




w 


W a 


Wl 


W2 


2 

Xmin 


ACDM 


265 +(u,i4 

U.ZOO_ Q 013 










424.911 


CPL 


U.ZOO_ Q24 


A - u '— 0.29 








423.432 


const Az 


fl O7Q + 0.1M 
0.033 


1 n+u.18 




Z.y:O_0.92 


-0.03 +U0ti * 


422.513 


const nAz 


0.2681^ 


±.UU_o 21 




1 7 i+o.yy 

L - 1 ± -1.41 


_ Q 41+0.43 


421.642 


free Az 


0.2761^ 


u - 9 ' -0.25 




L - oo -0.32 


-0 02 +u uti 

u - uz -3.39 


419.737 



* For the const Az model, the lower bound of W2 is only weakly constrained, so we do not list its value. 
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Figure 1. Marginalized probability contours at lcr and 2er CL in the f2 m — wo and wq — 
w a planes, for the CPL model. Constraints from SNLS3, BAO+CMB, BAO+CMB+iJ , and 
SNLS3+BAO+CMB+i/ arc all shown. Our results are consistent with the Fig. 6 of [22]. No 
evidence of deviations from the ACDM model is found. 




Figure 2. The evolution of w(z) and q{z) along with z for the CPL model. The result is consistent 
with the cosmological constant in the lcr CL. Moreover, q(z) is an increasing function of z, which is 
different from the result of [29]. 



4.1 The results of the CPL model 

The SNLS3 sample has been used to investigate the CPL model by Sullivan et. al in [22]. By 
performing a global fit, they obtained wq = —0.905^'^ and w a = — 0.984^}'^, which is con- 
sistent with the ACDM model. In this work, we obtained the similar results. In Fig. 1, we plot 
the la and 2a CL contours of the CPL model in the O m — wq and wq — w a planes, respectively. 
Constraints from SNLS3, BAO+CMB, BAO+CMB+# and SNLS3+BAO+CMB+# are 
shown in contours with differnt colors. The overlaps of these contours imply that these cos- 
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mological techniques lead to consistent cosmological implications in the CPL ansatz. We find 
that the result is consistent with the ACDM (wq = — 1 and w a = 0) at the la CL. This is 
different from the result of [29], where the authors found that the ACDM appears to be in 
tension with the observations (see Fig. 5 of [29]). 

To have a more direct view on the property of DE and the cosmic expansion, in Fig. 
2 we plot the evolution of the EOS w(z) and the deceleration parameter q(z) along with z. 
From the result of w(z), one can see that ACDM is within the la CL. Moreover, the evolution 
of q(z) shows that the current universe is under accelerating expansion. No evidence of the 
recent slowing down of the cosmic acceleration is found. This is also different from the result 
of [29]. The difference between the result of [29] and ours may arise from the difference of 
the data used in the two papers. 

4.2 The results of the three binned models 

In Fig. 3, we plot the evolution of w(z) for the three binned models. The best-fit as well 
as the la and 2a error bars are shown. From Fig. 3, it is clear that the error of w(z) in a 
certain redshift region depends on the number of SNIa in that region. The larger the number 
of SNIa, the smaller the error is. For example, the evolution of DE is most tightly constrained 
in the first redshift region (about a half of the SNIa in the SNLS3 sample lie in the redshift 
region z < 0.3), and is only weakly constrained in the last region at z > 1.0 (only 10 — 30 
SN lie in this region). This signature also manifests itself in the result of the free Az model, 
which has most tightly constrained w(z) in the second redshift bin. 

Except for slight difference in the values of best-fit and errors, the results of these three 
models are similar to each other, with the following common features: 

(i) The results of the three binned models are all consistent with the ACDM in the la CL. 

(ii) For the last redshift region, the upper bound of w(z) is constrained tow< 0.05, mainly 
due to the inclusion of the CMB data. The lower bound is only weakly constrained. 

(iii) At low redshifts, the result is tightly constrained and well consistent with the ACDM 
model. So the evolution of w(z) shows no evidence for the recent slowing down of the cosmic 
acceleration. 

To further confirm conclusion (iii), we also plot the evolution of q(z) for the three 
models at the low redshifts (z < 0.2) in Fig. 4. The results of the three binned models all 
indicate a universe with continuous accelerated expansion. This feature is in accordance with 
the result of the CPL model. 

In addition, we also plot the constraints on Q m for the these three binned models in 
Fig. 5. The figure shows that these three models lead to similar constraint on O m , with the 
best-fit value of 0.26-0.28, which is also consistent with the constraint on Q m in the ACDM 
model and the CPL model (see Table II). Therefore, we obtain consistent results from these 
four models. 

5 Conclusion 

In this work we probe the cosmic acceleration by using various cosmological observations, 
including the SNLS3 SNIa sample, the CMB data from the WMAP7 observations, the BAO 
data from the SDSS DR7, and the Hubble constant measurement from the HST. For models, 
we consider the CPL parametrization and three binned parameterizations based on different 
binning methods. We focus on the reconstruction of the DE EOS w(z) and the deceleration 
parameter q(z), to study whether the cosmic acceleration is slowing down, mentioned by 
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Figure 3. The evolution of w(z) along with z for the three binned models. From top to bottom, the 
results of the const Az, the const nAz, and the free Az models are shown. 



Shafiello et al. in [29]. Adopting four models in this paper, we obtain consistent results: the 
cosmological constant is still consistent with the data at la CL, and there is no evidence for 
the slowing down of the cosmic acceleration. Our result is different from the result of [29], 
which may arise from the difference of the data used in the two papers. Although our result 
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Figure 4. The evolution of q(z) along with redshift for the three binned models at the low redshift. 
From left to right, we show the results of the const Az, the const nAz, and the free Az models. We 
found no evidence for the slowing down of the cosmic acceleration. 




Figure 5. The distribution of the likelihood C oc e x I" 1 for the three binned models. From left to 
right, the results of the const Az, the const nAz, and the free Az models are shown. 



supports the standard ACDM paradigm, the uncertainties in the reconstruction of w(z) still 
allows space for other dynamical DE models. So the current observational data are still too 
limited to determine whether the cosmic acceleration is driven by a cosmological constant. 
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